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ABSTRACT 

Continuing advances in radar technology are driven both 
by new application regimes and technological innovations. 
A family of applications of increasing importance is detec- 
tion of a mobile target intruding into a protected area, and 
one type of radar system architecture potentially well suited 
to this type of application entails a network of cooperating 
radars. In this paper, we study optimal radar deployment for 
intrusion detection, with focus on network coverage. In con- 
trast to the disk-based sensing model in a traditional sensor 
network, the detection range of a bistatic radar depends 
on the locations of both the radar transmitter and radar 
receiver, and is characterized by Cassini ovals. Further- 
more, in a network with multiple radar transmitters and re- 
ceivers, since any pair of transmitter and receiver can po- 
tentially form a bistatic radar, the detection ranges of dif- 
ferent bistatic radars are coupled and the corresponding 
network coverage is intimately related to the locations 
of all transmitters and receivers, making the optimal de- 
ployment design highly non-trivial. Clearly, the detectabil- 
ity of an intruder depends on the highest SNR received by 
all possible bistatic radars. We focus on the worst-case in- 
trusion detectability, i.e., the minimum possible detectabil- 
ity along all possible intrusion paths. Although it is plau- 
sible to deploy radars on a shortest line segment across the 
field, it is not always optimal in general, which we illustrate 
via counter-examples. We then present a sufficient condition 
on the field geometry for the optimality of shortest line de- 
ployment to hold. Further, we quantify the local structure 
of detectability corresponding to a given deployment order 
and spacings of radar transmitters and receivers, building on 
which we characterize the optimal deployment to maximize 
the worst-case intrusion detectability. Our results show that 
the optimal deployment locations exhibit a balanced struc- 
ture. We also develop a polynomial-time approximation al- 
gorithm for characterizing the worse-case intrusion path for 
any given locations of radars under random deployment. 
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Figure 1: Bistatic radar SNR contours as Cassini 
ovals with foci at radar transmitter T and radar 
receiver R for distance products: c\ < ci < C3 < 
c 4 . 

1. INTRODUCTION 
1.1 Motivation 

An active radar system consists of a collection of 
transmitters and receivers in which transmitters emit 
RF signals and the receivers capture signals resulting 
from scattering of the transmitted signals reflected by 
objects of interest (targets). A monostotic radar is one 
in which a single transmitter and single receiver are col- 
located. A bistatic radar is also comprised of a single 
transmit-receive pair, but the transmitter and receiver 
are at different locations. Multistatic radar refers to a 
configuration that consists of multiple receivers at dif- 
ferent locations. Although physical-layer issues in radar 
have been extensively studied p], very limited atten- 
tion has been paid to radar network design 2 4. No- 
tably, [2] considered a network of monostatic radars, 
and [2111] studied a multistatic radar with one trans- 
mitter. We are thus motivated to consider a network of 
multiple radar transmitters and receivers where any pair 
of transmitter and receiver can form a bistatic radar. 
Unlike traditional passive sensors, one salient feature 
of radar sensing is that it actively emits signals to illu- 
minate the target. Further, radar technology provides 
superior penetration capability by using radio waves. 
Worth mentioning is that the cost of radars can be much 
higher than traditional sensors. 

The problem of detecting intrusion into a protected 



area, across a border or security perimeter, is of increas- 
ing interest for numerous applications, and has great po- 
tential in several regimes of growing importance, such 
as border monitoring, industrial and airport security, 
and drug interdiction. Worst-case coverage [5] is a com- 
mon criterion of efficacy in intrusion detection that has 
been applied extensively in the context of sensor net- 
works [6HS]. It is intimately related to the minimum 
possible detectability of an intruder (target) travers- 
ing a monitored field. In a bistatic radar network, the 
worst-case coverage is quite different and more diffi- 
cult to quantify than traditional sensor networks, be- 
cause 1) departing from the disk-based sensing range 
of a traditional sensor, the detection range of a bistatic 
radar depends on the locations of both radar transmit- 
ter and receiver and is characterized by the Cassini oval 
(in scenarios where intruder detectability is dominated 
by SNR). Formally, a Cassini oval is a locus of points 
for which the distances to two fixed points (foci) have 
a constant product (as illustrated in Figure [T]); 2) the 
sensing ranges of different bistatic radars are coupled 
with each other, since each transmitter (or receiver) can 
pair with other receivers (or transmitters) to form mul- 
tiple bistatic radars, indicating that its location would 
impact multiple bistatic radars. 

1.2 Summary of Main Contributions 

• We first rigorously quantify the intrusion detectabil- 
ity and network coverage attainable by a bistatic 
radar network, while taking into account the com- 
plication that each pair of transmitter and receiver 
can potentially form a bistatic radar. We then 
study the worst-case coverage under determinis- 
tic deployment, aiming to find optimal deployment 
locations of radar transmitters and receivers such 
that the worst-case intrusion detectability is max- 
imized. We present a sufficient condition on the 
field geometry under which it is optimal to de- 
ploy radars on a shortest line segment across the 
field, which we refer as a shortcut barrier. The 
line barrier appears intuitive but does not hold for 
arbitrary field geometry, for which we give counter- 
examples. Further, when radars are deployed on 
the shortcut barrier, the worst-case intrusion de- 
tectability turns out to be the vulnerability of the 
shortcut barrier, which is the minimum detectabil- 
ity of all points on it. 

• A main thrust of this study is devoted to character- 
izing the optimal deployment locations of radars 
on the shortcut barrier to minimize its vulnera- 
bility. Since the optimization problem is neither 
diffcrentiable nor convex, it is highly non-trivial 
to solve. We first quantify the local structure of 
detectability corresponding to a given deployment 



order and spacings along a shortcut barrier. Then, 
for a given deployment order, we establish the exis- 
tence and optimality of balanced deployment spac- 
ings required to attain the minimum vulnerability. 
Next, we derive sufficient conditions for an opti- 
mal deployment order, and characterize the cor- 
responding optimal deployment orders. Our find- 
ings reveal that the optimal deployment locations 
exhibit a balanced structure. Furthermore, under 
the optimal deployment, it suffices for a receiver 
to form at most two bistatic radars to guarantee 
the optimal coverage quality, under the condition 
that the number of receivers is no less than that 
of transmitters (which is often the case since radar 
transmitters are more costly). 

• We also study the worst-case coverage under ran- 
dom deployment, and quantify the worst-case in- 
trusion path for any given deployment of radars. 
In particular, by developing a novel 2-site Voronoi 
diagram with graph search techniques, we design 
an algorithm to find an approximate worst-case in- 
trusion detectability, where the approximation er- 
ror can be made arbitrarily small. The algorithm 
is shown to have polynomial-time complexity. 

We believe that the studies we initiated here on bistatic 
radar networks scratch only the tip of the iceberg. There 
are still many questions remaining open for the design 
of a networked radar system. 

The rest of this paper is organized as follows. Sec- 
tion [2] introduces the model of worst-case coverage for 
the bistatic radar network. In Section [3l we charac- 
terize the optimal deployment of radars for maximizing 
the worst-case detectability. Section [4] presents an ef- 
ficient algorithm for finding an approximate worst-case 
detectability given arbitrary locations of radars. Sec- 
tion [5] provides the evaluation results. Related works 
are discussed in Section [BJ The paper is concluded in 
Section [71 

2. SYSTEM MODEL 

Basic Setting. We consider a bistatic radar network 
deployed in a field of interest F . The field F is a 
bounded and connected region enclosed by four curves: 
a left boundary Fi, a right boundary F ri an entrance 
side, and a destination side (see Figure [2]). Note that 
we allow the boundary curves of the field to be ar- 
bitrary. An intruder (target) can traverse through F 
along any intrusion path P C F from a point on the 
entrance to a point on the destination. Specifically, the 
bistatic radar network consists of M radar transmitters 
Ti € T, i € M = {1,- ■ ■ ,M} and N radar receivers 
Rj e 11, j e Af = {1, • • • ,N}. We assume that each 
pair of transmitter and receiver can potentially form 
a bistatic radar. We further assume that orthogonal 



transmissions are used for interference avoidance. For 
convenience, we also use Tj or Rj to denote the location 
(point) of node Tj or Rj, respectively. 

Bistatic radar SNR and Cassini Oval. For a bistatic 
radar, the received SNR for a target is given by 

SNR = i|TX||WII 2 W 
where ||TX|| and ||PA| denote the transmitter-target 
and receiver-target distances, respectively, and K de- 
note bistatic radar constant which reflects physical-layer 
system characteristics, such as transmitting power, radar 
cross section, and transmitter/receiver antenna power 
gains. For ease of exposition, we assume the bistatic 
radar constant is the same, and our study on homoge- 
neous bistatic radars here will serve as a major step for 
studying the heterogenous case. Given a bistatic radar, 
the SNR contours are characterized by the Cassini ovals 
with foci at the transmitter and receiver. We use Ct,r(c) 
to denote the Cassini oval with foci at points T and R 
for constant distance product c. 

Intrusion Detectability and Network Coverage. Let 
||a&|| denote the (Euclidean) distance between two points 
a and b. It follows from (JTJ that the received SNR by 
a bistatic radar {Ti,Rj} from a target at some point 
p is determined by the distance product ||Tjp||||i?jp||. 
Then the detectability of a point p is quantified by the 
minimum distance product from p to all bistatic radars, 
denoted by 



lip) = min ||Tk>||||.R,-p|| 



(2) 



Observe from ([2]) that the detectability increases when 
I(p) reduces, and vice versa. The detectability of an 
intruder is captured by the detectability of its intrusion 
path P, which is the maximum detectability of all points 
on P, denoted by 



B(P) 



minPn). 

PGP 



(3) 



From the radar network's perspective, in the worst case, 
the intruder traverses through P along an intrusion path 
P* such that B(P*) is maximized. We refer this path 
and its detectability as the worst-case intrusion path 
and the worst-case intrusion detectability (WID), re- 
spectively. 

3. DETERMINISTIC DEPLOYMENT 

We consider both deterministic deployment and ran- 
dom deployment scenarios, in the same spirit as the 
arbitrary/random network models in the seminal work 
by Gupta and Kumar [10] , In the deterministic de- 
ployment case, our goal is to find optimal deployment 
locations of radar nodes (i.e., M transmitters and N 
receivers) in F such that WID is maximized, i.e., 

minimize B(P*) (4) 

Ti£F,Rj£F 
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Figure 2: (a) H is a shortcut barrier; (b) H is a 
shortest barrier but not a shortcut barrier since 
H' is shorter than H; the optimal deployment 
may not be on H. 

The above problem is highly non-trivial in general since 
the boundaries of field F can be arbitrary. In particular, 
we are interested in the line-based deployment scheme in 
which radars are deployed along a line which is simple 
to implement in practice. Further, we will show that 
it is optimal under some mild condition on the field 
geometry. To the best of our knowledge, the optimality 
of line-based deployment for worst-case coverage has not 
been studied in literature. 

3.1 Shortcut Barrier-based Deployment 

We need the notion of barrier to establish our re- 
sults on the line-based deployment. Given any c > 0, 
the detection range of a bistatic radar {Ti,Rj}, de- 
noted by AT it Rj (c), is the region enclosed by the Cassini 
oval CTi^jic) such that ||Pjp|| < c for all p e 

A-Ti^j (c). The detection range of the bistatic radar 
network, denoted by A(c), is the union of the detection 
ranges of all the bistatic radars. A barrier is a curve in 
the field F connecting the left boundary Fi and right 
boundary F r such that any intrusion path through F 
intersects with the barrier. The worst-case coverage is 
intimately related to the notion of barrier as given in 
the next property. The proof is based on an argument 
similar to the growing disks in [7] and is omitted here. 

Property 1. The worst-case intrusion detectability 
B(P*) is equal to the smallest value of c such that there 
exists a barrier in the network's detection range A(c). 

We define the vulnerability of a barrier U as the mini- 
mum detectability of all points in U, denoted by 



QiU) = max/fa). 

pG(7 



(5) 



Intuitively, we expect to attain the best coverage qual- 
ity by deploying radars along a shortest barrier which 
an intruder must pass through. We make the following 
assumption on the geometry of field F. 

Assumption 1. There exists a shortest line segment 
H connecting p and F r such that H C P. 

Although a shortest line segment connecting p and F r 
always exists, it is critical to assume that it lies in the 



field F (see Figure [2]) ■ We refer such an H as the short- 
cut barrier (SCB) and let h denote its length. In other 
words, h is the shortest distance between a point in Fi 
and a point in F r . Note that a shortcut barrier is always 
a shortest barrier that has the shortest length among all 
barriers, but the converse is not true (as illustrated in 
Figure [5]). Assumption 1 captures a large class of field 
geometry, e.g., any F of convex shape belongs to this 
class. 

The next result shows that the SCB-based deploy- 
ment is indeed optimal if the shortcut barrier exists. 
(All proofs of the theorems and lemmas in the sequel 
are relegated to Appendix) . 

Theorem 1. Under Assumption 1, it is optimal to 
deploy radars on the shortcut barrier H , in order to 
maximize the worse- case intrusion detectability. 

Although Theorem [1] appears intuitive, we caution 
that it is non-trivial since the proof hinges on the exis- 
tence of a shortcut barrier such that a line barrier with 
no greater vulnerability can be constructed from any 
arbitrary barrier. In other words, if the shortcut bar- 
rier does not exist, the optimal deployment may not be 
on a shortest barrier (see Figure 2(b)). 

Let Hi and H r denote the end points of H with 
Hi G Fi and H r G F r . Note that the worst-case in- 
trusion detectability under a line-based deployment is 
no greater than, but not necessarily equal to, the vul- 
nerability of the line segment. However, the equality 
holds for the SCB-based deployment as we show in the 
next result. 

Theorem 2. Under Assumption 1, if radars are de- 
ployed on the shortcut barrier H , the worst-case intru- 
sion detectability amounts to the vulnerability of H , i.e., 
B(P*) = Q(H). 

Combining Theorem [1] and [5J we conclude that the 
problem (|4]) is equivalent to finding the radar deploy- 
ment on the shortcut barrier H such that Q(H) is min- 
imized. 

3.2 Line-based Optimal Deployment Locations 

In this subsection, we study the optimal deployment 
locations of radars on the shortcut barrier H to mini- 
mize Q{H). Let ti = \\HiTi\\ and rj = \\HiRj\l Mathe- 
matically, our problem can be formulated as 



minimize 

ti ,Tj 



max mm \x 

0<x<hi£M,j&M 



(6) 



subject to < ti < h, i G M 
<r 3 <h,j G A" 

where mmi & M,jeAf \ x ~ U\\% — rj\ represents the de- 
tectability of a point p G H with \\Hip\\ = x. It can 
be easily checked that the objective function of prob- 
lem © is neither differentiable nor convex in general. 



Therefore, standard optimization methods can not be 
applied here. 

Consider the reformulation of problem (|6]) as follows. 
First, we treat Hi and H r as two virtual nodes and 
relax the constraint HiH r = h. Suppose all nodes in 
T and 1Z as well as Hi and H r are deployed on a line 
such that Hi and H r are the leftmost and rightmost 
nodes, respectively. Let S = (Hi,Si,--- ,Sj,H r ) de- 
note a deployment order ("order" for short) of all nodes, 
where J = M + N and (Si, ■■■ , Sj) is a permutation 
of the set T UK such that \\HiHi\\ = < \\HiSi\\ < 
■■■ < \\HiS.;\\ < \\HiH r \\. Without ambiguity, S also 
denote an order of locations of all nodes. Let Dg = 
(H-ff/Sill, • • • , ||Sjif r ||) denote the deployment spacings 
("spacings" for short) of a given deployment order of 
nodes S, which are the distances between all pairs of 
neighbor nodes in S. Similarly, let (Si, Si+i, • • • , Si+j) 
denote a deployment suborder ("suborder" for short), 
which is a deployment order of neighbor nodes in S, 
and T)(Si,s i+1 -,s i+i ) = (\\SiSi+i\\, ■ ■ ■ , \\S i+ j-iS i+ j\\) 
denote its deployment spacings. Clearly, any deploy- 
ment order S with any deployment spacings Ds repre- 
sent some deployment locations of T and K on a line 
segment with length \\HiH r \\. Likewise, any deployment 
locations of T and 1Z on the line barrier H can be repre- 
sented by some deployment order S with some deploy- 
ment spacings Ds under the constraint ||i?j.ff r || = h. 
To summarize, the problem ^ can be recast as 

Problem 1. Find an optimal deployment order S* 
with optimal deployment spacings Dg„ under the con- 
straint \\HiH r \\ = h such that Q(HiH r ) is minimized. 

In what follows, we outline the main steps to solve 
Problem 1. 

Stepl: We show a general and important structure of de- 
tectability on HiH r (LemmaQJ, which leads to the 
concept of balanced spacings. 

Step2: We introduce the concept of local suborder. Then 
we show the existence and characterize balanced 
spacings for a local suborder (Lemma [2]) , and es- 
tablish the optimality of balanced spacings for a 
local suborder (Lemma [3]). 

Step3: We rule out some orders from consideration, which 
leads to the concept of candidate order (Lemma|4|) . 
Then we show that a candidate order consists of 
decoupled local suborders, such that the results in 
Step 2 can be applied to show that balanced spac- 
ings exist and are optimal for a candidate order 
(Theorem O). 

Step4: Based on the optimal spacings obtained above, we 
derive sufficient conditions for an optimal order 
(Theorem 0|) . Then we characterize the optimal 
orders. 
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Figure 3: The detectability of a point attains 
local minimums at the end points and the mid- 
points of all pairs of neighbor radar nodes. The 
local vulnerable values are (a) unequal with ar- 
bitrary deployment spacings; and (b) equal with 
balanced deployment spacings. 

We start with two observations. First, swapping the 
locations of any pair of transmitters or any pair of re- 
ceivers results in an equivalent deployment. Second, 
transmitters and receivers are reciprocal to each other. 
Specifically, replacing all transmitters by receivers and 
replacing all receivers by transmitters results in an equiv- 
alent deployment. These observations will be used re- 
peatedly in the sequel. 

Let Y a b denote the midpoint between two points a 
and b. In the next lemma, we show a local structure for 
the vulnerability of HiH r . 

Lemma 1. Given any order S with any spacings Ds, 
I(p) attains local maximums on HH r at the end nodes 
Hi and H r , and at the midpoints of all pairs of neigh- 
bor nodes in (Si, ■ ■ ■ , Sj) (as depicted in Figure^ . In 
particular, 

arg max I(p) = Hi, arg max I(p) = H r 

pGH,Si P es.,H r 

arg max I(p) = Y SiSi+ , , Vz G {1, ■ • • , J — 1}. 

For convenience, we refer a local maximum point of I(p) 
and its value as a local vulnerable point and local vulner- 
able value, respectively. The local structure presented 
in Lemma Q] is essentially due to that the detectability 
of a point is dominated by the closet transmitter and 
receiver. Therefore, it suffices to characterize the local 
vulnerable values at the local vulnerable points to de- 
termine the vulnerability of HiH r . We say the spacings 
of an order of nodes are balanced if all the local vul- 
nerable values between or on the two end nodes of that 
order are equal (see Figure [3]). 

We define the pattern of an order of nodes as the or- 
der of node types (transmitter type T or receiver type 
R) of those nodes. We use T k or R k to denote k consec- 
utive T or R in a pattern, respectively. We define three 



kinds of local patterns Pi = (T,R k ,H r ), (R,T k ,H r ), 
(Hi,R k ,T), or (H h T k ,R) for k > 1, P 2 = (T, R k , T) 
or (R, T k , R) for k > 2, and P 3 = (T,R) or (R,T). 
We say a suborder Si is a local suborder if it has a 
local pattern, and it is associated with a local zone, de- 
noted by Zsi, which is the line segment between the 
two end nodes of that local suborder. For example, if 
Sj = (Ti, R u ■ ■ ■ , R k ,H r ), then Z Si = T\H r - The vul- 
nerability of a local zone Zs t is Q(Zs i ). We denote the 
length of Zsi by Ls ; . 

We observe that for any local suborder Si with any 
spacings, the closest transmitter and receiver from S 
for any local vulnerable point in Z^, t are nodes from Si . 
For example, if Sj = (Ti, R\), we can see that the clos- 
est transmitter and receiver for Yr lRl are Ti and R\, 
respectively; if Si = (T\,R\,--- ,Rk,H r ), the closest 
transmitter for any of Y TlRl , ■■ ■ , Y Rk _ lRk , H r , is Ti, 
and the closest receiver for any of Yr lRl , ■ ■ ■ , Y RklRk , 
H r , is from {R%, ■ ■ ■ , Rk}- Hence, the vulnerability of 
a local zone Zs t , i.e., Q(Zs t ), only depends on the lo- 
cations of nodes in Si, or in other words, depends on 
the spacings Ds 4 , and it is independent of the locations 
of nodes not in Si. Similar observations can be made 
for any local suborder. Using Lemma ([I]), in the next 
result, we characterize the balanced spacings of a local 
suborder. 

Lemma 2. For any c > 0, let eo(c) = 2y/c and ej(c) 
be the unique positive value of x such that (y~]?Z ei(c) + 
x/2)(x/2) — c for any j > 1. Given any c > and any 
local suborder Si, there exist unique spacings Ds ; such 
that T>Si is balanced with Q(Zsi) = c. Furthermore, it 
is given by, e.g., if Si = (Ti,Rx), then D Si = (e ); 
if Si = (Ti, Ri, ■■■ , R k , H r ), then 

Ds; = (eo, ei, • • • , e k -i,e k /2); 

if Si = (Ti, Ri, ■ ■ ■ , R k , T2) and k is even, then 

L>Si = (eo, ei, • • • , e k / 2 -i, e fc/2, e fc / 2 _i, • • • , e\, e ) 

or if k is odd, then 

Ds, = (eo, , e(fe_i)/2, e(fe_i)/ 2j •" , ei, e ) 

where the arguments "(c)"o/ej(c) are omitted for brevity. 
Similar results can be obtained for any other local sub- 
order. 

By definition, given c, the value of ei(c), i > can 
be found iteratively and it decreases as i increases (see 
Table 1). 

The following lemma shows that balanced spacings 
are optimal for a local suborder. 

Lemma 3. Given any local suborder Si, supposeDsi 
is balanced with Q(Zs t ) = c. IfS^ is the same local sub- 
order as Si and with spacings Dg' such that Q(Zs> ) < c, 
then Ls> < Ls { . Furthermore, if Ls> = Ls it then 
TV D s 



Table 1: Values of balanced spacings 



c 


eo(c) 


d(c) 


e 2 (c) 


e 3 (c) 


e 4 (c) 


1 


2.0000 


0.8284 


0.6357 


0.5359 


0.4721 


5 


4.4721 


1.8524 


1.4214 


1.1983 


1.0557 


10 


6.3246 


2.6197 


2.0102 


1.6947 


1.4930 


20 


8.9443 


3.7048 


2.8428 


2.3966 


2.1115 



Lemma [3] suggests that given a local suborder and a 
constraint on the local zone's vulnerability, the balanced 
spacings maximize the length of the local zone. 

We need the following lemma which reduces our search 
space significantly. 

Lemma 4. There exists an optimal order S not hav- 
ing a suborder with any of the following patterns: (T, T , 
R,R), (R, R,T,T), (T,T,R,H r ), (R,R,T,H r ), (Hi,T, 
R,R), and (H h R,T,T,). 

We say an order S is a candidate order if it does not have 
a suborder with any of the patterns given in Lemma 2J 
We observe that any given candidate order S consists 
of a series of local suborders Si,-- - , S m such that 1) 
each node in S is included in some S^; 2) the last node 
of Si is the first node of S,+i for all i = 1, - - • , m — 1. 
In particular, Si and S m have patterns Pi, and S; has 
pattern P 2 or P3 for i = 2, • • • , m — 1. For example, 

Si s 3 

S = (Hi, Hi, R2,Tx, R3, i?4, R5, T 2 , Re, T 3 , T4, H r ). 

" * ' " v ' 

S2 S4 

Accordingly, HiH r consists of a series of local zones Zs 1 , 
■ ■ ■ , Zs m , and Ds consists of a series of spacings Ds x , 
• • ■ , Ds m . To see why such a series of local suborders 
exists for any given candidate order S, we construct a 
super order S + from S by combining neighbor nodes 
of the same type in S into a super node. A node not 
combined in S is a simple node in S + . For the above 
example, the super order is given by 

S + = (Hi, Ri 2 , Ti, R 3j 4 i5 , T 2 , R e , T 3 + 4 , H r ). 

It is clear by our construction that two neighbor nodes 
in S + (excluding Hi and H r ) are of different types 
(transmitter type or receiver type). Since S is a candi- 
date order, it does not have a suborder with any pattern 
given in Lemma [4] Hence, it can be easily checked that 
two neighbor nodes in S + can not be both super nodes, 
and a super node can not be the third or third-to-last 
node in S + . Then we can see that the first three nodes 
and the last three nodes in S + represent two local sub- 
orders with patterns Pi; each super node and its two 
neighbor simple nodes in S + represent a local suborder 
with pattern P 2 ; two neighbor simple nodes not at the 
beginning or end of S + (excluding Hi and H r ) represent 
a local suborder with pattern P3. 



The decoupled structure of a candidate order allows 
us to apply results of local suborders. Consider a given 
candidate order S. We know that S consists of a se- 
ries of local suborders Si, • • ■ , S m . As a result, the set 
of local vulnerability values on HiH r is a union of dis- 
joint sets of local vulnerability values on all the local 
zones Zs i; ■■■ , Zs m . Therefore, since local vulnerabil- 
ity values on a local zone Z$, i only depend on the corre- 
sponding spacings Ds ; , we can see that Ds is balanced 
with Q(HiH r ) — c if and only if Ds; is balanced with 
Q(Zs i ) = c for all i = 1, • • • , m. Hence, by Lemma [51 
there exists unique Ds such that Ds is balanced with 

Q(RWr) = C. 

We observe that given a candidate order S, under the 
constraint that Ds is balanced with Q(HiH r ) = c, Ds 
varies as a function of c. It follows from Lemma [2] that 
each of ||if/S'i|j, • • • , ||5j_ff r || is an increasing function 
of c, and hence \\HH r \\ = WHS^ + ■■■ + \\SjH r \\ is 
also an increasing function of c. Furthermore, it can 
be easily verified that ||i?/if r || = when c = and 
— > 00 when c — ¥ 00. Thus, there exists a unique 
c > such that ||fijff r || = h. This implies that under 
the constraint ||7?;7? r || = h, there exist unique balanced 
Ds. We summary this result below as a corollary of 
Lemma O 

Corollary 1. Given any candidate order S, under 
the constraint \\HiH r \\ — h, there exist unique spacings 
Ds such that Ds is balanced. 

Applying Lemma [3l we obtain the following theorem. 

Theorem 3. Given any candidate order S, the unique 
balanced spacings Ds minimize Q(HiH r ) under the con- 
straint \\HiH r \\ = h. 

Theorem [3] suggests that the optimal spacings for any 
given candidate order is the unique balanced spacings 
referred in Corollary [TJ 

Since we have found the optimal spacings D s for any 
candidate order S, our next step is to find an opti- 
mal order S*. Suppose the number of transmitters is 
no greater than the number of receivers, i.e., M < N, 
WLOG. Since all transmitters are equivalent, we sup- 
pose that transmitters in S are indexed by their order 
in S such that < ||i?/Ti|| < • • • < ||H/T M || < \\HiH r \\. 
Define N = (n HlTl , n TlT2 , ■ ■ ■ , h Tm _ iTm , n TuHr ) where 
nTiT i+1 , TiHiTu and nr M H r denote the number of re- 
ceivers between Ti and Tj+i, between Hi and T\, and 
between Tm and H r , respectively, in S. Since all re- 
ceivers are equivalent, it suffices to determine nr ( r j+1! 
n^Ti, and Ut m h t for an optimal order. The next re- 
sult provides sufficient conditions for the optimality of 
an order. 



(a) Hi _Q □ Q Q Q □ Q_Q □ Q_ H T 



(b) Hi _Q □ OOP □ OOP □ Q_ H r 



(c) H, Q_Q □ Q O O □ P P P □ CL H r 



Figure 4: Optimal deployment locations of radar 
transmitters (square) and receivers (circle) for 

(a) M = 3 and N = 7 (b) M = 3 and N = 8 (c) 
M = 3 and TV = 9. 

Theorem 4. .A candidate order S is optimal if 

\n Ti T i+1 - n Tj T j+1 1 < 1, V« ^ j 

\ n TiT l+1 - 2nu l Ti \ < 1) l«TiT i+1 - 2nT M j? r | < l,Vi. 

Using the optimality conditions in Theorem^ we can 
find an optimal order S* described as follows. Let two 
integers q and r be the quotient and remainder of N/M. 
If q is even, there exists an S* such that 

r M-l-r 

N* = ,<7 + l,?~^79, |); 

if g is odd and r = 0, there exists an S* such that 



if q is odd and r > 1, there exists an S* such that 

r -! Af-r 

N = (— — ,q+l,--- ,q + l,q,--- ,q,— — ). 

It can be easily seen that any order obtained from the 
above optimal order by exchanging the values of 
and riT M H r , or the values of n.T 4 T i+ i an d «TjT j+ i for 
some i ^ j, also satisfies the optimality condition, and 
hence is optimal. We observe that the optimal order 
also exhibits a balanced structure. 

Given an optimal order S*, by Theorem [31 the opti- 
mal spacings D s » is the unique balanced spacings under 
the constraint ||iJ;i? r || = h. The optimal spacings can 
be found by a bisection search. In each search step, 
for a search criterion c, we can obtain by Lemma [2] the 
unique balanced spacings Ds* with Q(HiH r ) — c, and 
hence obtain \\HiH r \\. If \\HiH r \\ > h, c is decreased 
for the next step; if ||ff;H r || < h, c is increased for the 
next step. The search completes if ||-ff;iJ r || = h. Upon 
completion, the optimal spacings Dg„ is found and the 
optimal value of Q(HiH r ) is equal to the search crite- 
rion c. 



4. RANDOM DEPLOYMENT 

Next, we consider random deployment that can also 
be of great interest. In this scenario, given arbitrary lo- 
cations of radars resulted from random deployment, we 
aim to find the worst-case intrusion path P* and the 
worst-case intrusion detectability B{P*). Due to the 
complex geometry of bistatic radar SNR, it is difficult 
to find P* accurately. Therefore, we design an efficient 
algorithm for finding an intrusion path P' whose de- 
tectability is arbitrarily close to B(P*). Our algorithm 
partitions the field F into sub-regions based on a novel 
2-site Voronoi diagram, and then constructs a weighted 
graph G for the sub-regions to search for an approxi- 
mate worst-case intrusion path Pq . 

Given a set of points S, the Voronoi region Vs(s) of 
a point s € S is the set of points which are closer to s 
than to any other point in S, i.e., 

Vs(s)±{p\mm\\p S '\\ = \\ps\\}. (7) 

s' Go 

The Voronoi diagram for S is the collection of all Voronoi 
regions {Us(s),s G S}. To partition the field F, we 
first construct Voronoi diagrams for the set of trans- 
mitters T and the set of receivers 1Z, respectively, i.e., 
{Vr(Ti),Ti G T} and {V n (Rj),Rj G K}, within the 
boundaries of F. Then F can be partitioned into 2- 
site Voronoi regions, where a 2-site Voronoi region, de- 
noted by VfjiiTi, Rj), is the intersection of a transmit- 
ter Voronoi region and a receiver Voronoi region, i.e., 

VtM^, Rj) = V T {T t ) n Vn(Rj ) (8) 

for all T,eT and Rj e 11 such that V T {Ti) n Vn(Rj) ^ 
(see Figure [5]). Choose an arbitrary e > 0. For each 
2-site Voronoi region V-]-,n(Ti, Rj), we plot a sequence 
of Cassini ovals Ct { ,r (ke) for consecutive positive in- 
teger values of k such that Vq- ,n(Ti, Rj) is further di- 
vided by these Cassini ovals into sub-regions where each 
sub-region is bounded between two Cassini ovals with 
distance products being consecutive multiples of e. 

We construct a graph G — (V, E) where each vertex 
represents a sub-region, and two vertices are connected 
by an edge if they are adjacent sub-regions. Then we 
assign each vertex Vi a weight Wi, which is equal to 
the smaller distance product of the two Cassini ovals 
bounding this sub-region. We add two virtual vertices 
s and t to represent the entrance and destination of F, 
respectively. There exists an edge between vertex s and 
Vi for v-i £ V \ {s,t} if sub- region is adjacent to the 
entrance. Similarly, we add an edge between vertex t 
and Vi for Vi G V \ {s, t} if sub-region Vi is adjacent to 
the destination. 

It is clear that any intrusion path P is represented by 
an s — t path in G, denoted by Pa- We define the weight 
of an s - t path P G , denoted by W(P G ), as W(P G ) = 
mxa.u :Vie p a \ Wi. We aim to find an s — t path P' G with 
maximum weight. This can be obtained by a bisection 
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Figure 5: Overlapping Voronoi diagrams for 
radar transmitters and receivers. The polygon 
with indices {a, b, c, d, e] is a 2-site Voronoi region 
with respect to T\ and R\\ it is divided into sub- 
regions {v ll v 2l v Sl v i ,v b }. 

search between the smallest and largest vertex weights 
in G. In each step, breadth-first-search is used to check 
the existence of a s — t path using only vertices with 
weights larger than a search criterion c. If a path exists, 
c is increased to restrict the vertices considered in the 
next search step; otherwise, c is decreased to relax the 
constraint on the search. Upon completion, & s — t path 
with maximum weight is found. 

Suppose M — 0(n) and N — 0(n). The number of 
vertices in graph G, i.e., the number of sub- regions in 
field F, is 0(n 8 /e 2 + n 6 /e + n 4 ) and the reason is as fol- 
lows. We can treat the field F as a planar graph, where 
the vertices are the crossing points of Voronoi edges and 
Cassini ovals, and the edges are line segments or curves 
between two crossing points. The number of Voronoi 
edges and the number of Voronoi vertices are both 0(n). 
So the number of crossing points of Voronoi edges is 
0(n 2 ). Hence, the number of 2-site Voronoi vertices is 
0(n 2 ), and the number of 2-site Voronoi edges is 0(n 4 ). 
Since the field F is bounded, the number of Cassini 
ovals plotted for a 2-site Voronoi region is 0(l/e). So 
the total number of crossing points of Cassini ovals with 
2-site Voronoi edges is 0(n 4 /e). Hence, the number of 
vertices of sub-regions is 0(n 4 /e + n 2 ), and the number 
of edges of sub-regions is 0(n 8 /e 2 + n 6 /e + n 4 ). By Eu- 
ler's formula, the number of faces, which is the number 
of sub- regions, is 2 - 0(n 4 /e+n 2 ) + 0(n 8 /e 2 +?i 6 /e+n 4 ), 
equal to 0(n 8 /e 2 + n 6 /e + n 4 ). 

In the next result, part (a) follows from the definition 
of 2-site Voronoi diagram and our graph construction; 
part (b) follows from our result on the number of ver- 
tices in graph G. The proof is omitted due to space 
limitation. 

Theorem 5. (a) The maximum path weight obtained 
by our proposed approximation algorithm is within e to 



(a) H, Q DQ Q O n O Q CO Q Hr 



(b) Hi -Q □ Q Q Q □ Q Q Q □ Q_ H r 

Figure 6: Heuristic deployment locations of 
radar transmitters (square) and receivers (cir- 
cle) for M = 3 and N = 8 (a) HEU-1 (b) HEU-2. 

the worst-case detectability, i.e., 

B{P*)-e<W{P' G )<B(P*). (9) 

(b) Given e, our proposed approximation algorithm has 
polynomial-time complexity. 

5. PERFORMANCE EVALUATION 

In this section, we present some numerical results to 
illustrate the effectiveness of our proposed optimal radar 
deployment scheme. 

We examine the vulnerability of a line segment H un- 
der the line-based deployment. In particular, we com- 
pare the optimal deployment scheme (OPT) with two 
heuristic deployment schemes. The first heuristic (HEU- 
1) is to deploy radar transmitters (or receivers, respec- 
tively) with uniform spacings such that the maximum 
distance from a point on H to its closet transmitter 
(or receiver, respectively) is minimized (see Figure [5]). 
Specifically, HEU-1 results in 2||^Ti|| = \\TiT 2 \\ 
= ••• = \\T M -iTm\\ = 2\\T M H r \\ and 2||tf; J R 1 || = 
II-R1-R2II = • • ■ = ||i?jv-i-Riv|| = 2||i?Ar-ff r ||, where trans- 
mitters and receivers are indexed by their orders in S, 
respectively. It can be easily verified that this scheme is 
optimal if we treat radar transmitters (or receivers, re- 
spectively) as traditional sensors with disk-based sens- 
ing model. The second heuristic (HEU-2) is to deploy 
radar transmitters and receivers according to an op- 
timal order S* from OPT but with uniform spacings 
such that the maximum distance from a point on H 
to its closet transmitter or receiver is minimized (see 
Figure©. Specifically, HEU-2 results in 2|| J ff ; S' 1 || = 
1 1 5*52 1 = • • ■ = II 5£ f _ 1 5^ / 1| = 2||5^-.H r ||. 

Figure [71 Figure [S] and Figure M depict the vulnera- 
bility of H versus the number of radar receivers for 3, 
5, and 10 radar transmitters, respectively, where we set 
h = 100. We note that HEU-2 results in considerably 
lower vulnerability than HEU-1. This is because HEU-1 
deploys transmitters and receivers independently, while 
HEU-2 takes into account the joint design of transmit- 
ters and receivers. We observe that OPT further out- 
performs HEU-2 significantly. This suggests that opti- 
mal spacings is critical for improving the performance 
due to the complex geometry of bistatic radar SNR. 

Next, we evaluate the worst-case intrusion detectabil- 
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Figure 7: Vulnerability vs. number of radar re- 
ceivers for 3 radar transmitters. 
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Figure 9: Vulnerability vs. number of radar re- 
ceivers for 10 radar transmitters. 
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Figure 8: Vulnerability vs. number of radar re- 
ceivers for 5 radar transmitters. 
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Figure 10: Worst-case detectability vs. number 
of radar receivers for 3 radar transmitters. 



ity under random deployment using our proposed ap- 
proximate algorithm. We consider a square field F of 
size 100 x 100, where two opposite sides are entrance 
and destination, respectively, and the other two op- 
posite sides are left and right boundaries, respectively. 
We consider deploying transmitters and receivers in the 
field randomly (RAN) with uniform distribution. 

Figure HO] Figure [TT] and Figure [12] depict the worst- 
case intrusion detectability versus the number of radar 
receivers for 3, 5, and 10 radar transmitters, respec- 
tively, where we also plot the shortcut barrier-based op- 
timal deployment scheme (OPT). The results for RAN 
are averaged over 100 simulation runs. We observe that 
OPT performs significantly better than RAN. The rea- 
son is that deployment on a barrier is essentially much 
more efficient than deployment in the full field for worst- 
case coverage. 

6. RELATED WORK 

Worst-case coverage has been first studied in [5] for a 
traditional sensor network. Polynomial-time algorithms 
are devised to find maximum breach and maximum 
support paths between two locations. An efficient al- 
gorithm is proposed in [7J to solve the best-coverage 
problem raised in [5] . In [6] , efficient algorithms are de- 
veloped to find the minimum exposure path in sensor 
networks. Localized algorithms are designed in |llj to 
solve the minimum exposure path problem. A new cov- 
erage measure that captures both the best and worst- 
case coverage is studied in [5]. The deployment prob- 
lem to improve the maximal breach path is considered 
by [SUE]. 

Barrier coverage is an intimately related problem to 
worst-case coverage. The concept of weak and strong 



barrier coverage has been introduced in |13j , where crit- 
ical condition of weak barrier coverage is obtained for 
random deployment. The critical condition of strong 
barrier coverage is derived in [14] using percolation the- 
ory. An effective metric of barrier coverage quality is 
proposed in [T5]. [TB] studies constructing barrier by 
sensors with limited mobility after initial deployment. 
A novel full-view coverage model is proposed in |17] for 
constructing barrier in camera sensor networks. 

7. CONCLUSION 

Radar technology has great potential in many ap- 
plications, such as border monitoring, security surveil- 
lance. In this paper, we studied the worst-case coverage 
for a bistatic radar network consisting of multiple radar 
transmitters and radar receivers, where each pair of 
radar transmitter and receiver can form a bistatic radar. 
The problem of optimal radar deployment is highly non- 
trivial since 1) the detection range of a bistatic radar 
is characterized by Cassini oval which presents complex 
geometry; 2) the detection ranges of different bistatic 
radars are coupled and the network coverage is inti- 
mately related to the locations of all radar nodes. We 
present a general assumption on the field geometry un- 
der which it is optimal to deploy radars on a shortest 
line segment across the field, for maximizing the worst- 
case intrusion detectability. Further, we characterized 
the corresponding optimal deployment locations along 
this shortest line segment. Specifically, the optimal de- 
ployment locations exhibits a balanced structure. We 
also developed a polynomial-time approximation algo- 
rithm for characterizing the worst-case intrusion path 
for any given deployment of radars. 

To the best of our knowledge, the optimality of line- 
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Figure 11: Worst-case detectability vs. number 
of radar receivers for 5 radar transmitters. 
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Figure 12: Worst-case detectability vs. number 
of radar receivers for 10 radar transmitters. 

based deployment for the worst-case coverage, in par- 
ticular for bistatic radar networks, has not been studied 
before this work. Although the detectability model in- 
volves some idealized assumptions, we believe that this 
work will be of value in setting the foundations for net- 
worked radar systems. There are still many questions 
remaining open for the design of networked radar. 
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APPENDIX 

A. PROOF OF THEOREMS] 

Proof: Consider any given deployment locations {T, 7Z}. 
It suffices to show that there exist deployment locations 
{T',W} on H such that B(P*) > B'(P*'), where ' in- 
dicates that the deployment locations are {T',TZ'}. It 
follows from Property 1 that there exists a curve barrier 
U with U C A(B(P*)). Let U t and U r be the end points 
of curve U with Ui € Fi and U r <E F r . We will show 
that the deployment at {T',TZ'} forms a line barrier 
TWr with Q'tJWr) < Q{U) (see Figure EH). 

Let U' be the line passing through Ui and U r . Then 
we can find the projection of any Tj € T and Rj G 7Z 
onto line U' , denoted by T{ and R\ , such that T{T 1 i _L U' 

and RjRj -L U' . Consider any p 1 € UiU r . Since curve U 
is a path from Ui to U r , we can find a p £ U such that 
p 1 is the projection of p onto line V. Then it follows 
from the property of projection that \\Tip\\ > \\T(p'\\ and 
||-RjP|| > II^p'II f° r an h j- Hence, from ([2]) we have 
Hp) < I'ip'), and it follows that Q'(W7^) < Q{U) 
due to ([3]). Since H is the shortcut barrier, we have 
h < \\UiU r \\. Thus, we can translate the line-based de- 




Figure 13: An arbitrary radar deployment forms 
a curve barrier U (solid curve), while a line- 
based radar deployment can form a line barrier 
UiU r (dashed line) with no greater vulnerability. 



ployment {T',1Z'} onto the line where H lies such that 
Q'{H) < Q(U). Note that {V ,11'} may not be all on 
H. In this case, we can move each T[ ^ H or ■ ^ H to 
its closest end point of H such that I 1 (p) is not increased 
for all p £ H and hence Q'(H) < Q(U) still holds. 
Therefore, since H is a barrier and U C A(B(P*)), we 
have B'(P*') < Q'(H) < Q(U) < B{P*). □ 

B. PROOF OF THEOREM i 

Proof: Consider any given deployment locations {T, 1Z} 
on H. From the definition of barrier we have B(P*) < 
Q(H). Suppose B{P*) < Q(H). The proof is based 
on contradiction. It follows from Property 1 that there 
exists a barrier U with U C A(B(P*)) C A(Q{H)). 
We will show that barrier U must intersect with H (see 
Figure HH). 

Step 1: Consider a point p £ H with I(p) = Q(H). 
Suppose p ^ Hi and p ^ H r . Let H' be the line passing 
through p such that H' _L H. Consider any p' £ H' with 
p' ^ p. Since \\T t p'\\ > \\T lP \\ and H-R^'H > for 
all i,j, we have I(p') > I(p) and hence p 1 £ A(Q{H)). 
Thus H' intersects with A(Q{H)) at a unique point p. 

Step 2: Let Ai(Q(H)) and A r (Q(H)) denote the sub- 
sets of A(Q(H)) on the left and right side of H', respec- 
tively (see Figure [H|). Define O c (r) as an open disk 
centered at point c with radius r. We will show that 
MQ(H)) C Hl (h) and A r (Q{H)) C 0„ r (h)- 

Let Ohi (h) denote the circle centered at Hi with ra- 
dius h. Then let O l H[ (h) denote the subset of Ojj,(/i) 
on the left side of H', which is an arc. It can be verified 
that \\T z a\\ > \\Tip\\ and \\Rja\\ > \\Rjp\\ for all i,j and 

for all a £ Hi (h). Then since I(p) = Q(H), we have 

A(Q(H))nO l Hi (h) = 0. This implies that MQ(H)) C 
OHi(h). Similarly, we can show A r (Q(H)) C C>H T {h). 

Step 3: Let Ui and U r be the end points of U with 
Ui £ Fi and U r £ F r . Then we have Ui £ OH r (h) and 
U r O/f, (ft-), since otherwise \\UiH r \\ < h or ||fliJ7 r || < 
h, which contradicts that H is the shortcut barrier. Us- 
ing the result of Step 2, we have Ui £ Ai(Q(H)) and 
U r £ A r (Q(H)). Thus, curve U must intersect H'. 
From the result of Step 1 , U must pass through p. This 
implies that Q(U) > I(p) = Q(H), which contradicts 
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Figure 14: Any barrier with possibly lower vul- 
nerability than the shortcut barrier H must lie 
in the detection range A(Q(H)) (grey area) and 
pass through a line H' which has a unique inter- 
section p with A(Q(H)). 



that Q(U) < B{P*) < Q{H). 
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C. PROOF OF LEMMA [T] 

We start with an observations which is used repeatedly 
in this proof. Clearly, for any p £ H, I{p) only de- 
pends on the distances of p to its closest transmitter 
and receiver, say Ti and Rj, respectively. If p ^ TiRj, 
I(p) = \\Tip\\ \\Rjp\\ decreases as p moves towards Tj and 
Rj; if p £ T~R~, since \\T iP \\ + \\pRjW = \\T.RjW is a con- 
stant, it can be easily verified that I(p) = \\Tip\\\\pRj\\ 
increases as p moves towards Y^Rj, and hence attains 
maximum when p = Y^r ■ 

The main idea of this proof is to divide the line seg- 
ment between each pair of neighbor nodes into intervals 
such that all points on an interval have the same pair of 
closest transmitter and receiver, and then we examine 
the structure of detectability on each interval using the 
observation presented above. Due to space limitation, 
we only prove the case when Sk and Sk+i have different 
types (transmitter type and receiver type) for some k. 
The same idea is used to prove other cases. 

Suppose Sk and Sk+i have different types. Without 
loss of generality (WLOG), let Sk = Pi and Sk+i = Rj- 

1) Suppose the closest transmitter and receiver are Ti 
and Rj for any p £ TiRj. Since I(p) — \\Tip\\\\pRj\\ 
for p £ TiRj, I{p) increases as p moves towards Y^Rj- 

2) Suppose there exists some R m £ Hip such that 
Y RmR] £ TiY TiRj . Then I(p) = \\R m p\\\\ T lP \\ for p £ 
TiY RmRj , and 7(p) - \\T iP \\ \\pR 6 \\ for p £ Y RmRj Y TiRj . 
For each case of p, I(p) increases as p moves towards 
YT iRj - Similarly, if there exists some T n £ RjH r such 
that YfiTn € YpiRjRj, we can also show that I(p) in- 
creases as p moves towards Y TiRj for p £ YriRjRj- 

D. PROOF OF LEMMA H 

Due to space limitation, we only prove the case of Sj 
with pattern Pi, and the same idea is used to prove 
other cases. The main idea is to determine the spacing 



between two neighbor nodes successively according to 
the given order of nodes. 

Suppose Si = (Ti,i?i,-- - ,Rk,H r ) and Ds ; is bal- 
anced with Q(TiH r ) = c. It follows from I(Yr lRl ) = 
(||T 1 i? 1 ||/2) 2 = cthat||r 1 i? 1 || = 2^. Since I{Y Ri r 2 ) = 
\\TiY RlR2 \\ \\Y RlR2 R 2 \\ = {WTM+ HM2II/2) 
(||i?ii? 2 ||/2) = c, given \\TxRx\\ = we obtain a 

unique value of ||i£ii£ 2 ||. Similarly, in a recursive man- 
ner, given the values of ||Tii?i||, • • • , \\Ri-iRi\\, we ob- 
tain a unique value of such that I{Y RiRi+1 ) = 
c . . . until we obtain a unique value of ||i?A;-ff r || such 
that I(H r ) = c. Then we can see that ||Tii?i|| = eo(c), 
\\R1R2W =ei(c), \\Rk-xRkW =e fc _i(c), \\R k H r \\ = 
e/c(c)/2. Similar results can be shown for any Si with 
pattern Pi. 

E. PROOF OF LEMMA H 

Due to space limitation, we only prove the case of Si 
with pattern Pi and the same idea is used to prove 
other cases. Suppose Si — (Ti,i?i,--- ,R k ,H r ) and 
S' t = (Tl,R[,--- , R' k ,H' r ). Let I'(p) denote the de- 
tectability of p under S^. The proof is based on contra- 
diction. Suppose \\T[H' r \\ > \\TxH r \\. Since I(Y TlRl ) = 
(HT^ill/2) 2 = c > r(Y T{R[ ) = (\\T{R[\\/2)i, we have 
ll^i-Rxll < ll^i-Rill- Using this, we can show \\T[R' 2 \\ 
< ||7\i? 2 ||. Then following a similar argument, using 
\\T[R' 2 \\ < \\TxR2l we can show ||T^|| < ||Tii? 3 || ... 
until \\T[R' k \\ < \\TxR k \\. 

Since we have supposed HT^'if^l > ||Tiif r ||, we must 
havep^H = Hr^ll-IITXII > llTx^ll-HTAll = 
\\R k H r \\. It follows that I'{H' r ) = \\T[H' r \\\\R' k H' r \\ > 
\\TiH r \\\\R k H r \\ = c, which is a contradiction. Thus we 
conclude ||T;[if£|| < ||Tiif r ||. Further, it can be easily 
verified that if ||T{#;|| = ||Tiif r ||, we have \\T{R[\\ = 
||Ti.Ri||, \\R' k H' r \\ = H-Rfc^rll- Similar results can 
be shown for any Si with pattern P x . 

F. PROOF OF LEMMA H 

Suppose S = (Hi, • • • , Ti, T 2 , R%, R2, ■ ■ ■ , H r ) with any 
spacings. We can construct a new order of locations 
S' = (if/, • • • , Tl, i?i,T 2 ', R' 2 , ■ ■ ■ , H' r ) from S by swap- 
ping the locations of nodes T 2 and R\. Let I'(p) denote 
the detectability of p after the swapping. We observe 
that for any local vulnerable point p G ifjii, the clos- 
est transmitter is not T 2 , and for any local vulnera- 
ble point p € i2 2 if r , \\T2pW > \\T 2 p\\. Similar obser- 
vations can be made for Ri and R\. Thus we can 
see that I'(p) < I(p) for any local vulnerable point 
p € H1T1 or p e R 2 H r . We also observe that I(Yt 1 t 2 ) > 
(||TiT 2 ||/2) 2 = (HT^iH/2) 2 = I'{Y T , Rli ). Similarly, 
we can obtain I(Y R±R2 ) > I'(Yr> R >). Furthermore, 
we have I{Y T , Rl ) = (^RxW/Zf = (||i?i^||/2) 2 = 
/'(Y/j^). Thus we have shown that I'(p) < I(p) for 
any local vulnerable point p € HiH r , and hence, 
Q(HiH r ) is not increased after the swapping. Simi- 



lar results can be shown for an order with any pattern 
given in the claim by a swapping argument. 

G. PROOF OF THEOREMS] 

Suppose Ds is balanced with Q(HiH r ) = c. We have 
shown that S consists of a series of local suborders 
Si, ■ ■ • , S m . Suppose S' is the same order as S and with 
spacings D S ' such that Q'(H[H' r ) < c, where Q'(H[H' r ) 
denote the vulnerability of H[H' r under S'. Then S' 
also consists of a series of local suborders S'i , • • • , S' m 
such that Sj and S^- are the same local suborders for 
all j = 1, ••• ,m. It is clear that Q'(Zs>. ) < c, Vj. 
Then Lemma [3] implies that is'. < is j7 Vj. Since 

jyjLi Ls ., = h = J2T=i L S' S > we must have L s'j = Ls i' 
Vj. Hence, by Lemma [31 we have Ds' = Ds- 

H. PROOF SKETCH OF THEOREMS] 

Consider any given candidate order S. Clearly, S con- 
sists of suborders Si, ■ ■ • , Sm + i where Si = (if;, • • • ,Ti), 
• • ■ , Si = (7i-i, • • ■ ,Tj), ■ • • , Sjw+i = [Tm, ■ ■ ■ , H r ) . 
Suppose Ds is balanced with \\HiH r \\ = h and Q(HiH r ) 
= c. 

The proof has two steps. In step 1, we suppose the 
conditions given in the claim do not hold, and then we 
can find a new order S' from S with some spacings Ds' 
such that Q'(H',H r ) < c and H^'if^l > ||ifjif r ||, where 
Q'(H[H r ) denote the vulnerability of H[H' r under S'. 
This implies that S' is no worse than S. Repeating 
this argument, we can eventually find a new order S' 
satisfying the conditions in the claim and is no worse 
than S. In step 2, we show that all candidate orders 
satisfying the conditions in the claim are equivalcntly 
good. This implies that they are all optimal orders. 
Due to space limitation, we only prove step 1 for the 
case nTi_!Ti > nT j _ 1 T j + 2 and ny. > 2. The same 
idea is used to prove other cases in step 1 although more 
intricate arguments are needed for some cases. 

Suppose riTi^Ti > «t 3 _iT 3 + 2 and n^^Tj > 2. 
Clearly, Si and Sj must be local suborders. We can 
construct a new order S' from S by moving a receiver 
from between Ti_i and Ti in Si to between Tj-\ and Tj 
in S, . Then we set the spacings Ds' and Ds' balanced 
and Q(Z S >) = Q(Z S i ) = c, and we keep the spacings of 
any other suborder unchanged, i.e., Ds m = Ds^ for all 
m 7^ Therefore, we must have Q^Hj^H'.) = c. 

We observe from LemmaElthat we have is' — is = 

e («T J _ 1 T J +i)/2(c) if n Tj _ lT] is odd or e„ T ._ iT . /2 (c) if 
^Tj-iTj is even. Also, we have L Si -L s ^ = e ( „ T . _ iT ._i)/ 2 (c) 
if HTi-iXj is odd or e„ T lT ./ 2 (c) if n,T i _ 1 T i is even. Since 
eo(c) > ei(c) > e 2 (c) > • • ■ for any c, in any case, we 
have ||ff/fl-;||-||fl,ff r || = (i s; .-i Sj )-(i S( -is;) > 0. 
Recall that we have shown Q'(H[H' r ) — c. Then the de- 
sired result follows. 



